POSITIVE STABLE DENSITIES AND THE BELL-SHAPE 



THOMAS SIMON 

'^f^ • Abstract. We show that positive stable densities are bell-shaped, that is their n-th derivatives 

vanish exactly n times on (0, +00) and have an alternating sign sequence. This confirms the graphic 
predictions of Holt and Crow (1973) in the positive case. 
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1. Introduction 



Consider a smooth probability density denned on some open interval I C M and such that all its 
derivatives vanish at both ends of /. Rolle's theorem shows that the n— th derivative vanishes at 
least n times on /. The density is said to be bell-shaped if the n— th derivative vanishes exactly n 
times on / for every n > 1. For n = 1 this amounts to strict unimodality. For n = 2 this means, 
as for the familiar bell curve, that there is one inflection point on each side of the mode and that 
the second derivative is successively positive, negative, and positive. Rolle's theorem also entails 
that the zeroes of the successive derivatives of a bell-shaped density strictly interlace, and hence 
I/-} ■ that each derivative has an alternating sign sequence, starting positive. The sign sequence of the 

third derivative can be observed on a curve in comparing the aberrancy and the slope, as shown by 
£NJ \ Transon's formula [10], but contrary to the second derivative it requires an educated eye to guess 

o 

. what this sequence must be at first glance. The fourth derivative of a density has also a kinematic 

interpretation in terms of penosculating conies - see [10j and the references therein for a complete 
account, which however cannot be understood through the sole sequence of its signs. 

For a given density whose derivatives vanish at both ends of its definition interval, the usual way 
to show the bell-shape is to factorize its n-th derivative by a positive function and a polynomial 
function of degree n. This entails that this n-th derivative vanishes at most, and hence exactly, 
n times. Basic examples - whose details can be checked by the reader and serve as an exercise in 
calculus for sophomores - are the Gaussian, Gumbel and centered Student densities on R, and the 
inverse Gamma densities on (0, +00). When the density is not explicit, such factorizations are not 
appropriate anymore. Another criterion for the bell-shape is the ETP character of the associated 
additive convolution kernel - see Chapter 6.11.C in [7] - but the latter is rather difficult to check for 
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2 THOMAS SIMON 

non-explicit densities, and also quite stringent since it is e.g. fulfilled neither by the Student nor 
by the inverse Gamma densities. In this paper we are interested in the positive a— stable random 
variables X a and their densities f a , < a < 1. The latter are characterized by their Laplace 
transform and we choose the normalisation 

poo 

E[e- XXa ] = / e~ Xx f a (x)dx = e'^ , A > 0. 
Jo 

It is known - see e.g. the expansions (14.31) and (14.35) in [8j - that f a is real-analytic on (0, +oo) 
and that all its derivatives vanish at zero and at infinity. By Proposition 7.1.3 in [7] and (14.31) in 
[8], the kernel f a (% — y) is not TP2 and hence not ETP. The function /w 2 is explicit thanks to the 
identity in law 



X l/2 



d 1 



4Ti /2 

where here and throughout T a denotes the standard Gamma random variable of parameter a > 0, 
and is hence bell-shaped. The functions /1/3 and ^2/3 can be written down in terms of a modified 
Bessel resp. a confluent hyper geometric function - see e.g. (2.8.31) and (2.8.33) in [17J - but it 
does not seem to the author that the associated second order ODE's can provide any substantial 
information on their bell-shape. 

The strict unimodality of f a was first proved by Ibragimov and Chernin - see the proof of the 
Theorem in [6] - and then extended to all real stable densities by Sato and Yamazato - see Theorem 
1.4 in [S]. Since the latter are real-analytic, this actually follows from Yamazato's original theorem 
|15j on self-decomposable distributions. See also [12] for a short proof. The graphic simulations 
made by Holt and Crow [5], partly reproduced in [T7], display bell-shaped curves for all real stable 
densities, at least at the level of the second derivative. A consequence of Theorem 5.1. (ii) in [9] is 
that f a has only one inflection point on the left side of its mode, and more generally that each n-th 
derivative of f a vanishes only once before the first zero of the (n — l)-th derivative. The bell-shape 
property was claimed in [3] for all stable densities but a serious mistake on the TP character of 
the kernel f a (xy~ l ) - see Remark (f) in [12] for an explanation - invalidates this result, also at the 
unimodal level, except for a particular class of stable densities which will be discussed in Section 
3.3 and does not include the positive case. Apart from [5j [121 E] the bell-shape of stable densities 
seems to have escaped investigation, even in the most visual case n = 2. In this note, we show the 
full property in the positive case. 

Theorem. The densities f a are bell-shaped. 



Our main argument comes from Schoenberg's variation-diminishing property and the total pos- 
itivity of certain infinitely convoluted exponential kernels. The little known fact that X a admits 
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infinite exponential sums as additive factors had been noticed in |16j . It is used here in conjunction 
with a precise analysis of the other factor, which contains a certain exponential mixture with no 
atom at zero and can be chosen in order to have a weak bell-shape property of arbitrary order. 
In the case when 1/a is an integer, another proof of the bell-shape is obtained quickly by a mul- 
tiplicative factorization with inverse Gamma laws - see Section 3.1 below, but the argument does 
not extend to the other cases. With this method, the bell-shape of the density of log X\j n on R can 
also be established. In Section 3.2, we raise two natural conjectures on the bell-shape of positive 
self-decomposable densities. 

2. Proof of the theorem 

We first establish an additive factorization of X a via standard exponential variables. Recall that 
a positive random variable X ~ n is said to be an exponential mixture (ME) if its law has the form 

fj,(dx) = cSq (dx) + f(x)dx 

with c € [0, 1] and / a completely monotone function on (0, +oo). When c = we will use the 
notation X € ME*, and in this case Bernstein's theorem shows that the density / writes 

/•oo 

(1) f(x) = / 0e- 6 *n{d0) 

Jo 

over (0, +oo), where \i is some probability measure on (0, +oo). In other terms X = Exp (1) x X^ 1 
where here and throughout Exp (A) stands for the exponential law with parameter A, and X^ ~ fi. 
The following lemma is mostly due to Yamazato - see |16j pp. 601-602 - but we rephrase it and 
give a proof for the reader's convenience. 

Lemma (Yamazato) . One has the independent factorization 

(2) X a = Y a + Y ajTl , 

n>l 

where Y a>n ~ Exp((n7r/ sin(7ra)) 1//a ) for every n > 1 and Y a G ME*. 
Proof. We first write 

« r i i_ e -^ ) _ dx 

r(l - a) J x- ■ - jo V jo 

with the notation c a = sin(-7ra)/7r. This yields the further decomposition 

poo /'OO / roo 

X a = / (1 - e- Xx )l a {x)dx + / (l-e- Xx )( [(c a u) a }e- xu du)dx 
Jo Jo \Jo 

where [.] stands for the integer part and 

POD 

l a (x) = / {{c a u) a - [(c a u) a ])e- xu du 
Jo 



A« = T — — -J {l- e -*)—rz= (l-e- Xx ) {c a uYe-™du\dx 
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is a completely monotone function. By Theorem 51.10 in [8] there exists a random variable Y a 
belonging to the Bondesson class such that 

-logE[e~ AYa ] = / (1 - e- Xx )l a (x)dx, A > 0. 
Jo 

Because < {c a u) a — [(c a u) a ] < 1 for all u > and n~ 1 ((c Q ,u) Q — [{c a u) a ]) is integrable at zero, 
Steutel's criterion - see Theorem 51.12 in [8j - shows that Y a is actually ME. We next prove that 
Y a has no atom at zero, which will entail Y a E ME*. A change of variable and an integration by 
parts yield 

l a (x) = £. e -^ lu (^J\t a -[t a })dt^jdu 
for every x > 0, and it is a bit tedious but elementary to see that 

as u — > +oo. Theorem XIII. 5. 4 in [2] entails 

(3) l a (x) ~ g 

as x — > 0, so that Y a has infinite Levy measure and, by Theorem 27.4 in [8], no atom at zero. On 
the other hand, setting K Qjn = (nir/ sin(7ra)) 1/,Q! for every n > 1, one has 



roo poo 

/ [(c a u) a ] e~ xu du = J2 e~ xu du = 

•J JKrr.n. ~. ^ 1 



Hence, 



J (l-e- Xx )i^J [(c a u) a ]e' xu duj dx = ^"^^^ — dx 

= ^(log(A + K a , n ) ~ log(K a ,n)) 
n>l 

= -logE[exp[-A^Y a , n ]], 

n>l 

where the first equality comes from Fubini's theorem and the second from Frullani's. Putting 
everything together completes the proof. 

□ 

Remarks 1. (a) Because a € (0, 1), the sum on the right-hand side of ([2]) has finite expectation. 
This shows that X a and Y a have the same C p — integrability index which is a, excluded. The 
latter can also be seen from Theorem 25.3 in [8] and the fact that l a {x) ~ (a/T(l — a))x~^ a+1 ^ as 
x — > +oo. 

(b) Since the infinite exponential sum has a PF2 viz. log-concave density - see the proof of the 
theorem below for a stronger property - and since the density of f a is not PF2, one can interpret Y a 
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as the factor which breaks down the PF2 property for f a . Recall indeed that f a is log-concave until 
its second inflection point - see Theorem 1.3 (vii) in [9]. On the other hand Y a has a log-convex 
density by Holder's inequality. It is plausible that log/ a has only one inflection point, which can 
be readily checked for a = 1/2. 

We now introduce our notations concerning the sequence of signs of a real function defined on 
(0,+oo). We set 

{+ if x > ( ., . 

„ J + if n is even 

11 x = 1) and ± = < . . 

. r — 11 n is odd 

- if x < k 

for every x € [—00, +00], n £ N. For {u n } = {ex, . . . , e n } some finite sequence in {—,0, +}, we say 

that a continuous function / : (0, +00) — > R is of type E\ . . . e n if it has limits (finite or infinite) at 

zero and at infinity, vanishes on a finite set, and if the ordered sequence of its signs on [0, +00] is 

given by {ei, . . . , e n }. For brevity, we will write / ~ {u n } or / ~ e\ . . . e n to express this property. 

Observe that by the intermediate value theorem, the zero set of / ~ E\ . . . e n corresponds to a 

subsequence of period 2 in {e\, . . . ,£ n }- In particular, a function of type E\ ■ ■ ■ e n vanishes exactly 

(n/2 — 1) times on (0, +00) if n is even and either (n — l)/2 or (n — 3)/2 times on (0, +00) if n is 

odd. For example, the n-th derivative of a completely monotone function is of type ± n for every 

n > 0. The density function f a is of type 0+0 and, by the strict unimodality of X a , its derivative 

f' a is of type 0+0 — 0. Introducing the sequences 

{a n } = 0±°0± 1 0---± n and {b n } = ±°0± 1 ■ ■ ■ ± n 

for every n > 1, we can now state the central definition of this paper. 

Definition. For every n > 0, a smooth function f : (0, +00) — > R is said to be weakly bell-shaped 
of order n (WBS n ) iff® ~ {a t } Vi = 0...n and (-l)"+ 1 + i /W ~ {b n+1 } Mi > n + 1. 

The bell-shape property of a function / means that /W ~ {aj} Vi > 0, so that weakly bell-shaped 
functions are never bell-shaped. Notice also that WBS P H WBS,j = if p ^ q. The following result 
provides our key-argument. 

Proposition. Let X € ME* and < Ai < A2 < • • • < A n < . . . be some sequence. For every n > 1, 
the independent sum X + Exp(Ai) + • • • + Exp(A n ) has a WBS n _i density. 

Proof. Set fo for the density of X and recall from (TjQ) that it is completely monotone. In particular 
/o is real-analytic on (0, +00). We denote by f n the density of X + Exp(Ai) + • • • + Exp(A n ) for 
every n > 1. The latter are connected to one another by the formula 

(4) f n (x) = \ n e- X " x f e x "yf n ^(y)dy, 
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and hence solutions on (0, +00) to the linear ODE's 

(5) fn + Kfn = Kfn-l- 

An induction hinging on and © show that all densities f n are real-analytic on (0, +00), that 
all derivatives fn vanish at infinity, and that fn\o+) = for all n > i + 1 > 1. In particular one 
has f n ~ {00} f° r every n > 1. Using the complete monotonicity of fo for the initial step, another 
induction shows (— l) n+ Vn (0+) > for all i > n > 0, possibly with infinite values. This entails 
fn\o+) > for all n > 0, so that 0+ is an isolated zero of fn for all n > i + 1 > 1 (notice 
that this cannot be proved directly by analyticity because 0+ is also a singular zero). The weak 
bell-shape property will now be established through a more elaborate induction on n. 

We first show that f\ is WBSo- Since f\ ~ {a-o}, we know that /{ vanishes at least once on 
(0, +00). Set x\ = inf{x > 0, f[(x) = 0} > 0, the strict inequality coming from /{(0+) > 0. 
Differentiating (J5]) yields 

f'{{x) = Xtf^x) < 

for every x > such that f[(x) = 0. This shows that /{ vanishes only at x%, whence /{ ~ {^i}- 
Again, this entails that f" vanishes at least once on (0, +00). Set X2 = inf{x > 0, f'{{x) = 0} > 0, 
the strict inequality coming from /{'(0+) < 0. Differentiating ([5]) further entails 

f1'{x) = Ax/^x) > 

for every x > such that f'{{x) = 0. This shows that f" vanishes only at X2, whence — /{' ~ {&i}- 
Using the complete monotonicity of fo, the same argument yields (— ~ {61} for every 
i > 1, as required. 



We next show the induction step and suppose that f n is WBS n _i for some n > 1. We have 
already seen that f n +i ~ {^o} and we will first prove by induction on i that f n %x ~ { a i} f° r 
every 1 < i < n. We begin with f n+1 for the sake of clarity. Since f n +\ ~ { a o}i w e know that 
f'n+i vanishes at least once on (0, +00). Set x\ = inf{x > 0, f' n+1 {x) = 0} > 0, where the strict 
inequality comes from the fact that 0+ is an isolated zero of f n +i- Differentiating ([5]) as above 
entails 

Kf n ( x i) = fn+l(%l) < 0, 
where the inequality comes from the initial profile 0+0 of f n +i- Suppose that fn +1 (xi) = 0. Then 
by analyticity x\ is an isolated zero of so that 

hm — — — = -00. 
f n +l( x ) 
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By ([5]), this entails f' n {x) < as x ~\ x\. But the induction hypothesis yields f' n ~ {b\} if n = 1 and 
~ {ai} if n > 1, so that one cannot have f n (x{) = 0, a contradiction. Hence /4'+i( x i) < and if 

there exists X2 > x\ such that fn+i( x 2) = 0, then again the induction hypothesis entails f' n {x2) < 

so that fn+i( x 2) < 0, a contradiction. This shows that / ra+1 vanishes only once on (0, +oo) and, 

all in all, that f' n+ i ~ {<xi}. 

We now suppose f^+i ~ Wi} for some i < n and show /^!i ~ {a-i+i}- We know that /„!i 

vanishes at least i + 1 times on (0, +oo), and also from the above that has initial profile 0+0. 

Set {xj, j > 1} for the ordered sequence of its zeroes on (0, +oo). Again, one gets from ([5]) 

T J /n+1 0) T j /n+1 W 

for every j > 1. For clarity we consider two separate cases. 

(i) If i + 1 = n, the induction hypothesis entails fH +l) ~ {&„} = ±°0 - • ■ ± n 0. Set I* ++ {± fe 0} 
for the fc-th half-closed interval corresponding to each part ± fe of the graph of fn +l \ k = . . . n. 
The last equality in ([6]) with j = 1 and the initial profile 0+0 of /„!i entail that xi € i"* 1 with fei 
odd. If fn +1 \ x i) < 0, then Q with j = 1, 2 and the induction hypothesis show that X2 6 i^ 2 with 
^2 even. If fn +1 \x\) = 0, then X2 € l^p with &2 > k\. Iterating the procedure shows that there 
are at most, and hence exactly, n = i + 1 zeroes and that x^- € for every fc = 1 . . . n. All in all 
this shows that f^+P ~ { a i+i}j as required. 

(ii) If i + 1 < n, one has fn +l ^ ~ {oj+i} = 0±°0 - • • ± l+1 by the induction hypothesis. But 
adding a at 0+ in the profile does not change anything in the above analysis, and the fact that 
fn+i ~ { a «+i} follows exactly along the same lines as in (i). 

We finally show that (— l) n+1+i / n +i ~ {^n+i} for every i > n + 1. We only sketch the proof since 
the arguments are the same. Consider first the case i = n + 1. Since / n ™\ ~ {fln}> we know that 
/n+i^ vanishes at least n + 1 times on (0, +oo) and set {xk, k > 1} for the ordered sequence of its 
zeroes. Recall f!£^± 1 \o+) > and, from the induction hypothesis, that /i™ +1 ' > ~ ± 1 • • • ± n+1 0. 
Set again 1™ ^ {+ fc o} for the k-th half-closed interval corresponding to each part ± fc of the graph 
of fb n+1 \k = l...n + l. The same reasoning as above shows that necessarily x^ 6 7* for all k and 
hence that / n ™~^ vanishes exactly n + 1 times on (0, +oo), in other terms /„Ti ~ {^n+i}- It is 
then easy to show by induction on i that (— l) n+1+ */ n +i ~ {^«+i} f° r every i > n + 1, in using 
(_l)n+i+i^ i ( 0+ ) > o and the induction hypothesis {-l) n+l+i fH } ~ ± 1 0- ■ • ± n+1 0. 

□ 
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Remarks 2. (a) The complete monotonicity of /o is crucial in the above argument, initializing all 
induction steps. It seems difficult to extend the proposition to broader classes than ME. 

(b) When X itself has an exponential law and all parameters are different, the density of X + 
Exp(Ai) + • • • + Exp(A n ) is a linear combination of n + 1 distinct exponentials functions. Since 
the latter define a Chebyshev system - see [7j p. 24, the property WBS n is much easier to prove. 
We leave the details to the reader. 

(c) The method of examining the sign of the (n+l)-th derivative at the zero of the n-th derivative 
is used in [6] for n = 1 in order to show the strict unimodality of f a , after some explicit contour 
integration. The involved computations to extend this argument directly on the further derivatives 
of f a seem however extremely heavy, and the author does not believe that they could yield the full 
bell-shape. See [4J for an application of the contour integration of [6J to the asymptotic behaviour 
of these further derivatives. 

End of the proof. Fix n > 1 and set f a ^ n resp. g a ^ n for the density of 

n+2 oo 

X a + Y »,k resp. ^2 y Qjfc 

k=l k=n+3 

with the notations of the Lemma, so that one has 



/aO) = / fa,n( x ~ y)9a,n{y)dy. 
Jo 

By the Proposition, the function f a>n has a C n+1 extension on R and one obtains 

/■oo roo 

fi n+1) (x) = / f^ 1 \x-y)g a , n (y)dy = /& +1) (y)<? a ,n(x - y)dy. 



'0 ^0 

The Laplace transform of g a ^ n reads 

oo 00 



e Xx g a n(x)dx = 1 f ( — . 
o MJ fc lA 3 Vl + (sin(7ra)/fc7r)V«A 

and its reciprocal is of the Polya-Laguerre class £ * defined in [7] p. 336. By Theorem 7.3.2. (b) 
in [7J p. 345, this shows that the kernel g a ,n{x — y) is TPqo on R x R. Notice that g a ,n{x — y) is 
however not STPqo, because of the indicator function. Using the notation (3.1) p. 20 in [7] with 
I = (0, +oo), Theorem 3.1. (a) p. 21 in [7] and the Proposition entail 

S-(f< n+1) ) < S-(f a n +V) = n + l. 

Now since f a has isolated zeroes on (0, +oo) and vanishes at zero and infinity, Rolle's theorem 
yields 

S + (f a n) ) < S-(f a n +V)-l < n 
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where we have used the notation (3.2) p. 21 in [7J with / = (0, +oo). This shows that f& vanishes 
at most n times, and hence exactly n times, on (0, +oo). The proof is complete. 

□ 



3. Remarks and open questions 

3.1. The case when 1/a is an integer. In this situation the bell-shape follows easily from the 
independent factorization 

W *i/n = n~ n T~} n x ... x l£_ 1)/n , 

which was pointed out in [13] and also basically shown in [7] pp. 121-122. We discard the explicit 
case n = 2 and set g n resp. h n for the density of n~ n T~j n resp. Hr? x ... x rr^.. The 
multiplicative convolution yields 

fi/n( x ) = ( gn{xy~ l )h n (y) — 



and one can clearly differentiate under the integral. This entails 

fxU x ) = l°° gZHxy-^Kiy)^ = ^ K{xy~ l ) g ^{y)^- 1 d y 

for every x > 0, i > 1. The kernel /^(xy^ 1 ) is the composition of (n — 2) kernels that are STPqo on 
(0, +oo) x (0, +oo) by Theorem 2.1 p. 18 in [7j, and is hence itself STPqo on (0, +oo) x (0, +oo) by 
the Binet-Cauchy formula. With the above notation for S~ and S + , Theorem 3.1 (b) p. 21 in [7] 
and the aforementioned bell-shape property of g n yield 

S + (f[%) < S-(g$) = i 

for every i > 1 so that f\j n is bell-shaped, too. The above factorization (J7|) shows also easily that 
logXi/ n has a bell-shaped density on R, a property which we believe to be true for all a £ (0, 1). 
When a is rational, another multiplicative factorization of X a in terms of inverse Beta and inverse 
Gamma random variables was obtained in Lemma 2 of |11| . But the latter seems useless for the 
bell-shape of X a or logX Q , because inverse or log Beta laws have neither bell-shaped densities nor 
TPqo convolution kernels. 

3.2. Positive self-decomposable densities. A positive random variable X is self-decomposable 
when its Laplace transform reads 



E[e- xx ] = exp- 



70 A + / (l- e - Xu )^du 



A > 



u 

for some 70 > and k : (0, +00) — > R + non-increasing. For example, the positive stable law is 
self-decomposable with 70 = and k(x) = a/T(l — a)x a . It is known and easy to see - see Theorem 
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27.7 in [8] - that X has a density / when k ^ 0. The unimodality of / was shown by Wolfe 
In the following we will suppose that 70 = w.l.o.g. and set fk for the positive self-decomposable 
density associated with the non-degenerate spectral function k. 

When A;(0+) = +00, Theorem 28.4 (ii) in [8] and Lemma 2.5 in [9] show that fk is smooth and 
that all its derivatives vanish at zero and at infinity, and Theorem 1.4 in [9] entails that fk is strictly 
unimodal. When k{0+) < +00, it is also known - see Remark 28.6 in [8] - that fk is not smooth at 
zero and hence not bell-shaped. In view of the main result of the present paper, one might raise 
the 

Conjecture 1. The density fk is bell-shaped if and only if k(0+) = +00. 

Apart from the positive stable, examples of densities verifying the above conjecture are the 
inverse Gamma, or more generally the densities of T^~ a with t > 0, a > 1. Indeed, the latter are 
self-decomposable by the HCM criterion [I], and bell-shaped by a direct computation. Other self- 
decomposable examples are infinite exponential sums, whose bell-shape property can be proved 
exactly in the same manner as for f a . As pointed out in the introduction, when fc(0+) = +00 
Theorem 5.1. (ii) in [9] already shows that f^ vanishes only once before the first zero of 1 \ 
for every n > 2. The proof of this latter property relies on a repeated use of Steutel's equation - 
see e.g. Theorem 51.1 of [8]. Analyzing the further zeroes of fj n) with the same method seems 
however not quite obvious, because of the memory involved in this integro-differential equation. 

Let us now consider the case where k(0+) is finite. When fc(0+) < 1, Theorem 1.3 in |9] shows 
that fk is decreasing. When A;(0+) = 1, the situation is complicated because of the three different 
types I2T3 and I4 of [9]. In particular fk may have a non-trivial modal interval - see Lemma 53.2 
in [5J. When k(0+) > n + 1 for some n > 0, then Theorem 28.4 (i) in [8] and Lemma 2.5 in [9] 
show that fk is C n with derivatives vanishing at zero and at infinity, and again Theorem 1.4 in [9] 
entails that fk is strictly unimodal. The following is hence quite natural. 

Conjecture 2. // k(0+) > n + 1 for some n > 0, then ff® ~ {dj} for all i = . . .n. 

An example supporting this conjecture is the density f a ^ of — log f3 a j,, where (3 a ^ is the Beta 
random variable with parameters a > 0, b > 1. Indeed, one can show by a direct computation that 
f a ,b is self-decomposable with spectral function 

e ~ ax {\ — e~ bx ) 
= (l _ e -x) 

so that k a> b(0+) = b, and the explicit formula 

f h(x) - r(a + 6) ax( 
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entails that f a ^ satisfies the required properties. 

3.3. Two-sided stable densities. As mentioned in the introduction, stable densities are all visu- 
ally bell-shaped on Holt and Crow's graphics, and it is an interesting problem to show the property 
rigorously. It is known that two-sided stable densities are real-analytic on R, never vanish, and 
that all their derivatives tend to zero at infinity - see Remarks 14.8 and 28.8 in [8j. Hence, their 
n-th derivative vanishes at least n times on R by Rolle's theorem (a more complicated proof of 
this fact is also given in [3] pp. 234-235). The strict unimodality follows from Yamazato's theorem 
|15j by analyticity. The article [3] had claimed the full bell-shape property for all stable densi- 
ties with an argument relying on a certain integral representation through a Student-type kernel, 
and a multiplicative TP property. The integral representation is correct and basically amounts to 
Bochner's subordination and Zolotarev's duality - see Lemma 2 in [3j. However the multiplicative 
TP property claimed in Lemma 1 (iv) of [3j and crucial to bound the number of zeroes - see pp. 
236-237 in [3], is false in general. Indeed, Schoenberg's theorem shows - see pp. 121-122 and p. 
390 in |7J - that the kernel K a (x,y) = f a (e x ~ y ) is STPqo if and only if 1/a is the reciprocal of an 
integer. Actually, K a is not even TP2 for a > 1/2 - see the main result of |11| . In a forthcoming 
paper, we also show that K a is not TP n as soon as a > 1/n and 1/a is not the reciprocal of an 
integer. In [3] the skewed Cauchy case is also treated separately via Zolotarev's representation - 
see (2.8) therein. All in all the main result of [3] must hence be reduced to the 

Theorem (Gawronski). Two-sided a— stable densities are bell-shaped when a = 1,2 or when 1/a 
is an integer. 

Notice that this result does not extend directly to the one-sided case in spite of the affirmation 
made in (iv) p. 239 of [3], because the number of zeroes of a function sequence might clearly 
increase at the continuous limit. In order to cover all stable densities, the multiplicative method of 
[3] seems not appropriate because of the poor variation-diminishing properties of K a - see however 
|12j for a somehow related multiplicative point of view, which works at the unimodal level. On the 
other hand, Yamazato's additive factorization extends to all two-sided stable densities and delivers 
yet another proof of their unimodality - see again [16] pp. 600-601. Unfortunately the complete 
monotonicity of the factor Y a which is crucial in our argument to get the full bell-shape - see the 
above Remark 2 (a), is clearly lost in the two-sided case. 

Acknowledgements. Ce travail a beneficie d'une aide de l'Agence Nationale de la Recherche 
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